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5 ■ Abstract 

I ■ 
O ; 

I We perform a detailed analysis of the capabilities of the MINOS, ICARUS 

CIhI and OPERA experiments to measure neutrino oscillation parameters at the 

p) ' atmospheric scale with their data taken separately and in combination. MI- 

1^ '• NOS will determine Amfa and sin^ 26*23 to within 10% at the 99% C.L. with 
10 kton-years of data. While no one experiment will determine sin^ 2^13 with 

. ^ [ much precision, if its value lies in the combined sensitivity region of the three 

^ I experiments, it will be possible to place a lower bound of 0(0.01) at the 



d ' 95% C.L. on this parameter by combining the data from the three experi- 

ments. The same bound can be placed with a combination of MINOS and 
ICARUS data alone. 

PACS numbers: 14.60.Pq, 14.60.Lm 
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I. INTRODUCTION 



The hypothesis that neutrino oscillations are responsible for the atmospheric and 
solar pj-p^ neutrino anomalies is becoming increasingly favored as more data is accumu- 
lated and new experiments start running. In particular, the disappearance of atmospheric 
z/^ is currently being tested by K2K, the first long baseline experiment, with confirming 
indications ||13[|. The best-fit oscillation parameters relevant to the atmospheric and solar 
scales from the analysis of the data from the above experiments are (Amgg, sin^ 26'23) = 
(2.9 X 10-3 eV^ 0.98) and (Am^^, sin^ 2^12) = (4.9 x 10-5eV^0.79) ]15[, respectively. 
Similar results were obtained in other global analyses of solar data; see Ref. The remain- 
ing mixing angle, ^13, is constrained by the CHOOZ experiment to be small, sin^ 2^13 < 0.1, 
for values of Am32 relevant to atmospheric neutrino oscillations at the 95% CL. |[T7[ . 

If oscillations in the three neutrino framework are assumed to be the correct solution to 
the solar and atmospheric neutrino problems (although the participation of a sterile neutrino 
in the oscillation dynamics, motivated by the LSND |jl8| data, is convincingly consistent with 



all existing data |T9|), the mission of future experiments will be to precisely measure the 
relevant mixing parameters. The question of how well this task can be accomplished by 
the upcoming long baseline neutrino facilities is a rather important one since it will guide 
studies and set the goals for future neutrino experiments. 

The KamLAND reactor experiment |2^ will take important steps to solve the solar 
neutrino problem by determining if the large mixing angle MSW solution is correct, and by 
pinning down Am^i and sin^ 26*12 to accuracies of O(10)% if it is pi . 



Low statistics in the K2K experiment will limit its power to measure the oscillation 
parameters relevant to the atmospheric neutrino problem. Three other long baseline ac- 
celerator neutrino experiments will address the atmospheric neutrino deficit with much 



higher statistics: the Main Injector Neutrino Oscillation Search (MINOS) experiment [p^-|24 



from FNAL to the Soudan Mine, the Imaging Cosmic And Rare Underground Signals 
(ICARUS) experiment [25,26| and the Oscillation Project with Emulsion-t Racking Appa- 
ratus (OPERA) P7| , p5| from CERN to Gran Sasso. All are expected to start taking data 
around 2005, which will be long before the advent of superbeams or neutrino facto- 
ries for more references see Ref. |^|. These future long baseline experiments and 
KamLAND will therefore provide independent and accurate information about neutrino os- 
cillation parameters at the solar and atmospheric scales with the important advantage that 
their results will not rely on assumptions about the solar model or the atmospheric neutrino 
fluxes. 

In this work we explore the precision with which Am|2, sin^ 26*23 and sin^2^i3 can be 
measured by MINOS, ICARUS and OPERA in the framework of three neutrino oscillations. 
We include all oscillation channels accessible to these experiments and take into account 
the available experimental information on backgrounds and efficiencies. All three experi- 
ments will operate with a high i/^ flux and measure z/^ disappearance. Because of the energy 
dependence of the oscillation probability, the most sensitive measure of the oscillation pa- 
rameters is the energy spectra of the neutrinos reaching the detector. Observations of an 
energy dependent distortion of the spectrum can uniquely determine the oscillation param- 
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eters. Charged current events alone can provide information about the original z/^ energy 
spectrum. 

Although these experiments can also observe z/e, and to a lesser extent v^- appearance, 
the expected statistics in these channels are unlikely to be high enough to allow the energy 
spectra of these channels to be studied extensively. One can nevertheless extract important 
information from the total expected rates in these channels. 

The paper is organized as follows. In Sec. we very briefly review the three generation 
oscillation formalism. In Sees. [T] and |IV| we describe the MINOS, ICARUS, and OPERA 
experiments defining the observables relevant to our analyses and how they are calculated. 
In Sec. we outline the procedure used to analyse the data simulated for MINOS, ICARUS 
and OPERA. We also perform a global analysis of the data from the three experiments. In 
Sec. ^we discuss our results on the sensitivity to and precision with which Am32, sin^2^23 
and sin^2^i3 are expected to be determined by each experiment, and from the combined 



data of all three experiments. We present our conclusions in Sec. [VII 



II. 3 V OSCILLATIONS 

In the three generation framework there are three mixing angles and one CP phase that 
determine the amplitude for transition from flavor state to another. For the experiments 
under study, the CP phase has a negligible effect on the oscillation amplitudes and we set 
it to zero in what follows. The transformation between the neutrino mass eigenstates i^i, 1^2 
and z/3 and neutrino interaction eigenstates z/g, and v-r is given by 



C12C13 ■S12C13 Si3 

-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 
S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13 




(2.1) 



using the standard parameterization for the mixing matrix U [^. Here Cjj and Sij denote 
the cosine and the sine of the mixing angle 6ij. 

The propagation of neutrinos through matter |Q is described by the evolution equation 
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dr 
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Am^2 




(2.2) 



where Am?- = mf — and E, 



^^^^ is the neutrino energy. Here, Ae{r) = 2^2 Gpn^ 
^ eV^ Ye p (g/ cm^) Ey (GeV) is the amplitude for z/g — e forward scattering in mat- 



{r)E,, 



Here, Ae{ 

1.52x10"" 

ter with Yf, denoting the electron fraction and p the matter density. For the experiments 
considered in this work, the neutrino path only traverses the Earth's crust which has an al- 
most constant density of p ~ 3 g/cw? with Y^. ~ 0.5, giving A^, ~ 0.23 x 10~^ eV^ E^, (GeV). 
Although analytical expressions for the oscillation probabilities in matter of constant density 
exist 0], we will solve Eq. ( p.2|) numerically taking into account the dependence of density 



on depth using the density profile from the Preliminary Reference Earth Model . 
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Since the baseline of the experiments under consideration is 730 km, and lAmgj^l <^ 
lAm^gl, the contribution to the neutrino dynamics from the solar scale is small and the 
oscillation probabilities can be described in terms of just three parameters: Amgg, sin^ 2^23, 
and sin^2^i3. For a constant matter density one can approximate the probability that a 
muon neutrino is converted to an electron neutrino as []33| , ^ 



sm' 



2^]^ sin' 



23 sm^ 



V AE 



S 



(2.3) 



where 



S 



\ 



A, 



Amia 



— cos 26'i3 



+ sin^ 2^13 



(2.4) 



and 



sin^ 2^^ 



sin^ 29 
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(a^- COS 2^,3) +sin2 2^ 



(2.5) 
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For L ~ 730 km, matter effects make a negligible contribution to the probability for conver- 
sion into a tau neutrino and the vacuum expression serves as a good approximation, 



sin 2^23 cos ^13 sin 



' AE 



(2.6) 



From Eq. ( |2.3| ) one can see that the conversion probability for —> y^, is approximately 
proportional to sin^ 2^13 and is therefore small unless there is a resonant enhancement at 
cos 2^13 = Af./ Am\2- The conversion probability is, on the other hand, proportional 

to sin^ 26^23, and is therefore large. The survival probability, P^.^^^^ = 1 — Pu^^ur ~ Pu^^ue 
depends almost entirely on sin^ 2^23 and Am|2- Thus the parameters sin^ 26'23 and 
can be studied via z/^ disappearance or i>r appearance, while Ue appearance is necessary 
to probe sin^26'i3. The parameters sin^ 26*23 and Am32 should be readily accessible in all 
three experiments given the high event rates in the —>■ channel with z/,- detection also 
possible in both ICARUS and OPERA. Because of the expected smallness of sin^ 26*13 and 
the z/e background in the beam, z/e appearance will be difficult to study, and the accessibility 
to sin^ 26*13 depends strongly on the detector. 



III. MINOS EXPERIMENT 



The MINOS experiment [^] is designed to detect neutrinos from the Fermilab NuMI 
beam. The source of the neutrino beam is the decay of pions and kaons produced by 
collisions of 120 GeV protons (originating from the Fermilab Main Injector) with a nuclear 
target. A total of 3.7 x 10^° protons on target are expected per year. The beam will be 
almost exclusively z/^ with a (z/g + z/g) contamination of about 1%. A low, medium, or high 
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energy neutrino beam can be realized at MINOS by adjusting a focusing horn at the source. 
The resulting beam energies are peaked at approximately 3, 7 and 15 GeV, respectively. 
For MINOS to operate as a v^- appearance experiment either the high or medium energy 
beam would be required, since the r production threshold is 3.1 GeV. The experiment will 
most likely start with the low energy configuration, since this configuration will maximize 
its sensitivity to lower values of Arrigg, that are favored by the latest Super-Kamiokande and 
K2K results [jl4[. We do not include the effects of a hadronic hose or beamplug as these 



options are unlikely to be realized. 

There will be two iron/scintillator detectors associated with the MINOS beam. The 
1 kton near detector, which will detect neutrinos before oscillations occur, will be located 
on site at Fermilab. The 5.4 kton far detector will be located 732 km away from the source 
in the Soudan mine in Minnesota. It is expected that the experiment will run for two to 
three years starting in 2005. 

MINOS will independently measure the rates and the energy spectra for muonless (O/x) 
and single- muon (1/i) events, which are related to the neutral current (nc) and charged cur- 
rent (cc) reactions, respectively [|36|. The O/i and 1/i event rates can be used to measure the 



ratio of neutral current events to charged current events which is an important consistency 
check on neutrino oscillations and could be useful for determining if a sterile neutrino is in- 
volved in the disappearance channel. However, it does not enhance the precision with which 
the oscillation parameters can be extracted. Since our goal is to determine the precision 
with which the oscillation parameters can be found in a three neutrino framework, it is a 
sound assumption that z/^ disappearance is a consequence of transitions to Vr- We do not 
consider 0/i events in what follows. 

We assume the low energy beam configuration to study the MINOS sensitivity to Am|2, 
sin^ 26*23 and sin^ 26*13. We focus on z/^ — >• disappearance to calculate the 1/x-event en- 
ergy spectrum of charged current events and z/^ z/g appearance to calculate the total 
integrated rate of charged current z/g interactions. 

The 1/i-event energy spectrum is divided into 23 bins of variable width Ai?,, according 



to the method described in Ref. [0. The 1/x-event sample at MINOS will consist mainly 
of cc-events, with a small contribution coming from the misidentification of nc-events as 
cc-events. Thus, the content of the i^^ bin, dNl^/dE, is given by 



^(Am^2,sin2 2^23, sin^ 2^13) = ^ KiE)[yl{E)P,^^,^{E)e'^\E) 

+ a:^{E)v-%E)] dE . (3.1) 

Here, (j)^^ is the neutrino fiux at the MINOS far detector. We use the energy spectrum 
for z/^ cc-events in the MINOS far detector in Ref. |3^. The number of active targets is 
A = M 10^ Na Up Hy, where M is the detector mass in kton, 10^ A^^ is the number of nucleons 
per kton {Na = Avogadro's number), Uy is the number of years of data taking and Up is 
the number of protons on target per year. The z/^ charged current cross section, o"^^, and 



the neutral current cross section, af,'^, are provided in Ref. |Q and Ref. respectively. 



Here, e'^'^ is the probability that a given charged current Up, event will be correctly labelled 
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(~ 80% at the spectral peak), and ■i]'^^ is the probabihty of misidentifying a given neutral 
current event as a l/x event (~ 10% at the spectral peak). These factors were estimated by 
the MINOS collaboration using a detailed Monte Carlo simulation of the detector with the 
low energy beam configuration The Standard Model expectation at the far detector is 
about 400 single-muon events per kton-year. 

The MINOS experiment will reduce the systematic errors associated with the charged 
current energy spectrum by comparing the spectrum measured at the far detector with that 
measured at the near detector from the same beam. The uncertainties in the far/near ratio, 
mainly due to the theoretical uncertainties in the secondary production of tt and K in the 
NuMI target, range from 1-4% in the peak of the low energy beam up to 5-10% in the high 
energy tail |3^. We have neglected these errors in our study. 

The z/e rate is defined as 

Re{Aml^, sin^ 2^23, sin^ 2^13) = + S^^"^ + B^^^^^ + 5f + S^'^^) , (3.2) 
where the signal, Se, is 

S,{Aml„sm^ 2023, sm' 26^3) = aJ <P,^{E) P,^^,^{E)a:i{E) e^dE . (3.3) 



Here, = 0.28 is the signal efficiency and a^^ is the Uf, charged current cross section pS 



The electron background event contribution from the decay of tau leptons in the detector is 
S(^^)(A"^32,sin2 2^23, sin^ 2^13) = ^/ Ki^) P.^^uAEKliE) e^^dE , (3.4) 



where ej^^^'' = 0.066 is the reduction efficiency and aj;^ is the z/,- charged current cross 
section ||3^. The background from e//i misidentification is given by 

Bi'^'^^Aml,, sin' 2023, sm^2e,s) = Aj 0.^(E) P.^^,^(E)<^(E)r/('^'^)rfE, (3.5) 

where rj^i^f^^ = 0.001 is the misidentification probability. The background coming from the 
Ue beam contamination is 

E^^-)(Am^2, sin^ 2^23, sin^ 2^13) = aJ <P,^{E) P.^^,^(E)<(E) ef ^'"^ciE, (3.6) 

where gj^'^'^™) = 0.079 is the reduction efficiency. The final background included, B^'^\ from 
the decay of neutral pions created by neutral current interactions, is 

i?(-)(Am^2, sin2 2^23, sin^ 2^13) = ^ / Ki^) K%E) ^^rfi?, (3.7) 

where e^^^ = 0.016 is the reduction efficiency. The efficiencies used in calculating both the 
signal and the backgrounds for MINOS events can be found in Ref. ||41|| . 
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IV. CERN-GRAN SASSO EXPERIMENTS 



A new facility under construction at CERN will direct a z/^ beam 732 km to the Gran 
Sasso Laboratory in Italy where it will be intercepted by two massive detectors, ICARUS 
and OPERA. The number of protons on target at the CERN SPS source is expected to be 4.5 
X 10^^ per year, and the z/^ beam will have an average energy of 17 GeV. The fractions z/g/ 1'^, 
^fil^fi and z/,-/z/^ in the beam are expected to be as low as 0.8%, 2% and 10~^, respectively. 



The ICARUS p5[ detector will use liquid argon for its detection medium, and is expected 
to have an initial effective volume of 3 ktons with a 10 year running time. We conservatively 
assume an exposure of 20 kton-years because more data than this does not significantly 
improve either the reach of the experiment or the precision with which the oscillation pa- 
rameters are determinable. OPERA ||2^, another detector at the Gran Sasso consists of 
lead plates interspaced with emulsion sheets, and is expected to have an effective volume of 
2 kton and a running time of 5 years. 

We investigate the capabilities of ICARUS and OPERA as disappearance, Ve appear- 
ance, and Vr appearance experiments. We calculate the full energy spectrum for the charged 
current events, and rates for the charged current Ve and Vr events. Low statistics in the 
Ve and Vj. appearance channels makes a study of their full energy spectra unfeasable. The 
z/^ scattering energy spectra in ICARUS and OPERA consists of 16 bins of width AE = 2.5 
GeV. The number of events in the i^^ bin, dNJ[^/dE, is calculated as 

-^(Am^^, sin^ 2^23, sin^ 2^13) = "ttt / Ki^) P.,^u,{E) e{E) dE, (4.1) 

ail/ lAh/ JEi-1.25 GeV 

where (p^^ is the flux of z/^ arriving at the Gran Sasso laboratory and o"^^ is the charged 
current cross section for z/^-nucleon scattering, both of which are available in Ref. . 



Again, A is the number of active targets in the detector, and we assume a constant efficiency 
e = 0.98 (0.94) for ICARUS (OPERA) according to Refs. and @, respectively. If 

no oscillations occur, 2180 (2070) charged current events per kton-year are expected at 
the ICARUS (OPERA) detectors. 

The Uj. production rate at ICARUS, Rr, is 

Rr{Aml^, sin^ 2^23, sin^ 2^13) = + Bi^''^ + E^^--), (4.2) 

where the signal, S^-, is defined as 

5',(Am22,sin2 2^23, sin^ 2^13) = ^/ 0.,(i?) P.,^.,(i?)<(i?) dE , (4.3) 



and we use an overall efficiency e,- = 0.06 P3| , ^ . Here we have assumed that the r iden- 
tification at ICARUS will be made via the r — e leptonic decay only [^,^,|^. The 



overall efficiency therefore includes the branching ratio which is ~ 18% and the selection 
cuts efficiency. 

Since tau events will be detected via an electron final state at ICARUS, z/g backgrounds 
in the beam must be considered for this experiment. One source is z/^ — ^ z/g oscillations, and 
is given by 
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5(^^)(Am^2,sm2 2^23, 8111^2^13) = aJ P.^^,,(i5;)<(E) ej^^) dE . (4.4) 

Another important background contribution comes from the intrinsic z/g beam contamina- 
tion, and is given by 

i?(^^^-)(Am^2, sin^ 2^23, sin^ 2^13) = aJ P.,^.,(E)<(E) e^^"^^ dE . (4.5) 



The reduction efficiencies are e^^'^'* = 0.15 and '^""^ = 0.009 Hil. 

For the OPERA experiment, we have assumed r detection via its one-prong decay into 
leptons (/) and hadrons (h) EH]. Backgrounds can be safely neglected, and the integrated 



rate of charged current u^- interactions at OPERA is given by Eq. (|jD with = 0.087. 



This efficiency value includes Y.i=i,h Br(r i) as well as the selection cuts efficiency . 

The expression for the z/g production rate at ICARUS and OPERA, i?e, is identical to 
the one presented in Eq. ( p.2[ ) with the appropriate fluxes, cross sections, and efficiencies. 
For ICARUS, ee = 0.75, e^^ = 0.14, r^^^^) = 0.002, e^"'"'^^ = 0.19 and ej"'^ = 0.01 so as 
to emulate the effect of the selection cuts of Ref p6|. The corresponding parameters for 



OPERA are eg = 0.7, e^^^^ = 0.13 (which reproduces the electron detection efficiency of 
Ref. [|7|), 7]^'"'^ = 0.002, ej^^"^""^ = 0.19 and ef^^'^ = 0.016. 



V. DATA SIMULATION AND ANALYSIS 

In this section we describe the procedures used to determine the sensitivity and contour 
plots for the three experiments and their combination. 



A. DEFINITIONS 

1. Uf, 

For this oscillation channel we use the information from the charged current energy 
spectrum of each experiment. All three experiments will observe large numbers of events 
of this type, so that the number of events in each bin obeys Gaussian statistics. The 
function is then defined as 

xliAmI,, sin^ 2^23, sin^ 2^13) = E i^t' ' (^f ^ " (5-1) 

where A^*^ denotes the theoretical prediction for the number of events in bin i (for a given 
set of oscillation parameters), calculated according to Eq. ( p.l| ) for MINOS and Eq. ( [4.1|) 
for ICARUS and OPERA. N°^^ denotes the "observed" number of events in bin i. The 
simulation of the data giving N°^^ will be described in Section C. The number of bins is = 
23, 16 and 16 for /3 = MINOS, ICARUS and OPERA, respectively. For MINOS, the error 
matrix is defined as afj = SijN°^^ + (0.02YN°^^Nj^^, where the off-diagonal elements refiect 
a 2% overall flux uncertainty [^]. For ICARUS and OPERA we only consider statistical 
errors, yielding af, = 6ijN°^^. 
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2. 

Gaussian statistics will also be realized in this channel. Even in the case of no oscillations 
the background contributes a sufficient number of events to Rf. to warrant a Gaussian 
function, 

CDobs pth \ 2 
^ ^ M . (5.2) 

Here, the theoretical electron rate B}^ is given by Eq. ( ^.21 ) with the appropriate experimental 
parameters (fluxes, efficiencies, etc) for the experiment under consideration. For MINOS, we 
define the error as a\ = Rf'^ + (5™^ x Rf^^)'^ assuming a global systematic error 5™™ = 0.1. 
For ICARUS, we define the error as = Rf^^ + (^^yg^ x i^(t'eam)^2^ where only the beam 
contamination has been considered as a source of significant uncertainty with 6^^^^ = 0.05. 

For OPERA, we use = ^ + (5°^^! x 5(^'''"'^)^ + (^sTst-ii x ^^'^)^ where Bj^^^"^'^ and 
are the OPERA backgrounds produced by the Ue beam contamination and the misidentified 
neutral current interactions, respectively. Their associated systematics errors are given by 
Cyst-i = 0-1 Cyst-ii = 0-2, as suggested in Ref. [0. We do not include a systematic 
error for the misidentified neutral current interactions at ICARUS because it has not been 
estimated. We neglect the errors arising from other components of the background because 
they have not been determined on account of their expected insignificance. 



3. Ur 



The number of detected events in this channel is likely to be small for both ICARUS 
and OPERA and so we use a function consistent with a Poisson distribution. 



X^(Am22,sin2 2^23,sin2 2^ 



13 J 



(Rf' - Rf) + Rf' In 



R. 



obs ^ 



(5.3) 



Here, the theoretical rate i?*^ is calculated via Eq. ( [4.3|) , using the appropriate experi- 
mental parameters for ICARUS and OPERA. When the number of Ur event rates is greater 
than 5, we use a standard Gaussian function given by 



X^(Am22,sin2 2^23, sin^ 2^13) 



Dobs D' 

— -ft, 



(5.4) 



The errors considered in this case are purely statistical. 



B. REGIONS OF SENSITIVITY 



For the z/^ v^. channel, we calculate the sensitivity regions for MINOS, ICARUS 

32 J 



and OPERA in the (sin^ 26'i3, t\m\^ plane. The procedure used is to set R?^^ 
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i?e(Am|2, sin^ 2^^23) 0) with sin^ 2^23 = 0.85Q or 1, and minimize the function given in 
Eq. ( |5.2| ) with respect to the parameters sin^2^i3 and Am^g- Curves are then drawn at 
= "~ Xmin = 4.61 corresponding to the 90% C.L.. The sensitivities are calculated 
for each experiment separately as well as for their combined potential by using the same 

procedure with = Xl(mmos) + Xe{icarus) + Xe(opera)- 

We determine the sensitivity to sin^ 26'23 and Am^g in the — > channel for MINOS 
and in the z/^ z/^ channel for ICARUS and OPERA. For MINOS, we calculate N^^^ = 
(iA'^{^/(i-E'(Am|2, 0, sin^ 2^13) with sin^ 2^13 set either to or 0.1, its minimum or maximum 
allowed value. Equation ( p.l| ) is then minimized with respect to sin^ 2^23 and Am|2, and 
the contour of constant A^^ at the 90% C.L. is found. For ICARUS and OPERA, we set 
i?°bs ^ i?^(Am22,0,sin2 2^i3) in Eq. (|^)/(|^), again with sin^ 2^13 = or 0.1, minimize 
with respect to sin^ 2^23 and Am32 and trace contours of constant Ax^. For the combined 
sensitivity of all three experiments, the same procedure is repeated with a global defined 

'^'^ X X^(minos) ~^ -^r(icarus) ~^ -^r(opera)" 



C. PREDICTING THE ALLOWED REGIONS 



To predict how accurately MINOS, ICARUS and OPERA will determine the neutrino 
oscillation parameters at the atmospheric scale, we simulate data sets for the three experi- 
ments assuming oscillations occur with 



Aml^ = 3 X 10"^ eV^ , Amj^ = 5 x 10~^ eV^ , 



sin^ 2^23 = 1 , sin^ 2^12 = 0.8 , sin^ 2^13 = 0.05. (5.5) 

All parameters other than ^13 are close to the best-fit values from the current data. The value 
of ^13 is chosen to lie within (or close to) the sensitivity regions of all three experiments under 



consideration while obeying the CHOOZ limit |T7]. We have assumed a normal hierarchy 
(Am^2 > 0). Note that since matter effects are small at 730 km, the results for an inverted 
hierarchy will not be substantially different from that for a normal hierarchy. 

Each simulated data point is generated by randomly choosing a point from a Gaussian 
or Poisson distribution associated with the number of events expected theoretically. This 
procedure is followed for the 23 bins of the MINOS z/^ energy spectrum as well as the 16 
bins of the ICARUS and OPERA spectra, giving the N"^"" of Eq. This is also 

done for the charged current rates of all three experiments, yielding the Rf'^ in Eq. { f).2\j , 
and the Ut charged current rates at ICARUS and OPERA, giving the of Eq. ( |5.3D or 
( ^.4|) . In the case of the Ue and Ur rates where backgrounds are taken into consideration, the 
backgrounds and signals are calculated and simulated separately before being combined to 
give the rates used in the analyses. 



^This value is approximately the smallest value of sin^ 2^23 allowed at the 99% C. L. by a combined 
analysis of SuperKamiokande and K2K data |l^. We conservatively assume that ^23 < vr/4. If the 
converse were true, the sensitivity to sin^ 2^13 is greater than for sin^ 2^23 = 1- 
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The allowed regions are calculated for the individual experiments by combining the 
information available from each channel and defining total functions as follows: 

Xmmos(^"^32' sin^ 2^23, slu^ 2^13) = + Xe . 

xLrus(^"^32. sin^ 2^23, slu^ 2^13) = xl + xl + xl , 

X^pera(A^32, siu^ 2^23, siu^ 2^13) = xl + xl + xl ■ (5.6) 

The x^ functions are then minimized by varying Am32, sin^ 2^23 and sin^ 2^13 with /S.m^^ = 
5 X 10~^ eV^ and sin^ 2^12 = 0.8. Contours of Ax^ = 6.25 and 11.36 are made corresponding 
to the 90% and 99% C.L., respectively. These confidence level regions are projected on to 
two-dimensional subspaces of the three parameter space. 

To investigate whether combining data from the three experiments improves the deter- 
mination of the oscillation parameters, we define a global function as 

Xminos+icarus+opera(^'^32 5 2^23) siu 2^^i3) = Xminos ~^ Xicarus ~^ Xopera- (5-7) 

This x^ function is minimized and contours are made as described above. 

VI. RESULTS 
A. MINOS 

In Fig. we present the sensitivity (at the 90% CL.) of MINOS in the z/g appear- 
ance and z/^ disappearance channels for an exposure of 10 kton-years. In the Uf, appear- 
ance channel (left panel of Fig. |l]), we see a small decrease in the sensitivity of MINOS as 
sin^ 26'23 is decreased within its allowed range. MINOS is sensitive to sin^ 26'i3 > 0.05 for 
Am|2 ~ 3 X 10~^ eV^ and maximal z/^ — Ur mixing. For the z/^ disappearance channel, one 
can see that the region of sensitivity is not affected by the variation of sin^ 2^13 within its 
allowed range. MINOS is sensitive to Am|2 ^ 5 x 10~^ eV^ at maximal sin^26'23. The 
sensitivity of the neutral current to charged current event ratio to the leading oscillation 
parameters has been analysed in Ref. That sensitivity is comparable to the sensitivity 
of the i>n disappearance channel in Fig. |] Direct evidence for transitions to Ur is a very 
important aspect of the MINOS experiment, which can be accomplished by comparing the 
O/i and 1/i event rates. 

For the oscillation parameters of Eq. ( |5.5| ), the signal (background) is comprised of 2770 
(0) cc events and 15 (41) events. Note that we have included the contribution from 
misidentified nc-events in the z/^ event sample. In Fig. ^we show the allowed regions obtained 
by simulating data for the MINOS detector. The Ax^ contours are shown at the 90 and 
99% C.L.. Because of high statistics and an optimal L/E combination, MINOS should be 
able to pin down Amgg and sin^ 2^23 quite precisely, as is shown most clearly by the plot 
in the (sin^ 26'23, ^^32) plane. For 10 kton-years of exposure, we find that MINOS should 
be able to measure these parameters to within 10% at the 99% C.L.. However, MINOS will 
not make a determination of sin^ 2^13. Low statistics in the z/g channel are expected because 
of the smallness of ^13, and relatively large backgrounds of z/g are expected to be present in 
the beam. 
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B. ICARUS 



The sensitivity of the ICARUS experiment to the z/g and Vj. oscillation 

channels are shown in the left and right-hand panels of Fig. ^ respectively. 20 kton-years 
of exposure is assumed, and the contours are drawn at the 90 and 99% C.L.. For the 
z/g appearance channel, we find that ICARUS can access values of sin^2^i3 > 0.03 for 
Am|2 ~ 3 X 10^^ eV^ and maximal — z/^ mixing. From the appearance channel, we 
find that ICARUS is sensitive to Amgg > 1 — 2 x 10^^ eV^ at maximal mixing. 

The theoretical input (Eq. (U)) yields 41960 (0), 34 (380) and 27 (11) z/^, z/g, Vr signal 
(background) events, respectively. Figure | shows 90 and 99% C.L. contours obtained by 
simulating data for all three channels that can be studied by the ICARUS experiment by 
observing variations in = Xe + + Xr- From the three panels displayed in the figure, it 
can be seen that ICARUS can determine Amgg and sin^ 2^23 to within 30% at the 99% C.L. 
but will not provide a meaningful determination of sin^ 2^13. 

C. OPERA 

Figure |^ shows the sensitivity of the OPERA experiment to oscillation parameters at 
the atmospheric scale via the z/^ — > Ve and z/^ — > Vr channels, respectively. This plot is made 
for 10 kton-years of exposure, and represents the 90% C.L.. Using appearance, we find 
that OPERA will be sensitive to values of sin^ 26^13 > 0.2 for Am|2 ~ 3 x 10"^ eV^ and 
maximal z/^ — mixing, but this region is ruled out by CHOOZ. Exploiting Vj. appearance, 
OPERA is shown to be sensitive to Amgg ^ 1 x 10~^ eV^ at maximal mixing independent 
of the value of sin^ 2^13. Although OPERA has no considerable background for r events as 
well as a higher detection efficiency than ICARUS, the latter has a greater sensitivity in this 
channel because of its larger volume and consequently higher event rate. 

The number of z/^, z/g, z/^ signal (background) events expected for the oscillation param- 
eters used are 19870 (0), 14 (230) and 16 (0), respectively. In Fig. || we present contours of 
Ax^ at the 90 and 99% C.L. for OPERA. 

D. Global Analysis 

In Fig. 1^ we present the 90% C.L sensitivity regions after combining the data simulated 
for the MINOS, ICARUS and OPERA experiments. In the (sin^ 26^13, Am^g) plane, the 
combined sensitivity does not show much improvement over ICARUS alone, resulting in a 
sensitivity to sin^26'i3 > 0.02 for Am32 ~ 3 x 10~^ eV^ and maximal z/^ — Vr niixing. The 
sensitivity of the experiments in the (sin^ 2^23? ^^32) plane depends only slightly on the 
value of sin^ 2^13 between and 0.1. The combined sensitivity of the three experiments is 
Am|2 > 4 — 5 X 10~^ eV^ at maximal mixing which is the same as that of MINOS. 

In Fig. 1^ we show the allowed regions obtained from the combined analysis at the 90%, 
95% and 99% C.L.. The precision with which the combination of the experiments can de- 
termine Am32 and sin^ 2^23 does not improve from the analysis with MINOS alone, since 
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changes in Xminos dominates over any changes in the overall function. These parameters 
can therefore be determined to within 10% once again. However, the ability of these ex- 
periments to constrain sin^ 2^13 becomes feasible by combining them into a single analysis. 
For the data we simulated, we obtain sin^ 26*13 ^ 0.01 at the 95% C.L.. Similar results for 
sin^ 26*13 are obtained by setting = Xminos + xfcarus; since only MINOS and ICARUS are 
sensitive to values of sin^ 26*13 < 0.1 for Aml2 = 3 x 10~^ eV^. 



VII. CONCLUSIONS 

We have studied the potential of the future long baseline experiments ICARUS, MINOS, 
and OPERA, to ascertain the neutrino oscillation parameters at the atmospheric scale as 
separate experiments, and in combination. By simulating data at Am^g = 3 x 10~^ eV^, 
Am|i = 5 X 10"^ eV^, sin^ 2^23 = 1, sin^ 2^12 = 0.8 and sin^ 2^13 = 0.05 for all three ex- 
periments, and calculating the allowed regions given by this data, we have estimated the 
precision with which these experiments can determine the parameters Am|2, sin^2^23; and 
sin^ 26*13. We took into consideration z/^ disappearance and z/g appearance for all three exper- 
iments and appearance for ICARUS and OPERA. The precision with which the leading 
oscillation parameters can be determined by the combined data of the three experiments is 
illustrated in Fig. |. 

The range of values of Am^g and sin^ 26*23 allowed by the Super-Kamiokande atmospheric 
neutrino data will be significantly narrowed by the MINOS experiment. With an optimal 
(L/E) ratio for studying these parameters, MINOS should pin them down to within 10% 
at the 99% C.L. with 10 kton-years of data. See Fig. |. ICARUS and OPERA, though 
not as sensitive to these parameters as MINOS, will provide an important check on the 
neutrino oscillation hypothesis by observing tau events via the z/^ — * Ur channel. The value 
of sin^ 26*13 will be difficult for these experiments to determine because of low event rates 
and large backgrounds in the z/g channel. We have shown that this parameter will remain 
unbounded when the experiments are analysed separately. Combining the data from the 
MINOS and ICARUS experiments may allow a lower bound on sin^ 26*13 of (9(0.01) to be 
placed at the 95% C.L. if sin^ 2^13 lies within the combined sensitivity of the experiments. 
Establishing a lower bound on sin^ 2^13 would eliminate models that predict smaller values 
of sin^2^i3. If sin^ 2^13 > 0.01 and the CP phase is sufficiently large, CP violation in the 
lepton sector could be studied at superbeam facilities 

In the next decade, with data from the three experiments considered (and a little luck), 
we could have good knowledge of the parameters that mediate atmospheric neutrino oscilla- 
tions. Along with KamLAND's possible determination of Amfg and sin^ 26*12 using reactor 



neutrinos all the elements of the mixing matrix other than the CP phase could be 
known. 



13 



ACKNOWLEDGMENTS 



Wc thank A. Bettini for a critical reading of the manuscript and A. Erwin, M. Goodman 
and D. Harris for discussions. We are grateful to M. Messier for providing us with the MINOS 
energy spectrum for the low energy beam and to D. A. Petyt for information on MINOS 
efficiencies. D. M. and B. P. W. thank Universidade de Sao Paulo for its hospitality during 
the initial stages of this work. This work was supported by Fundagao de Amparo a Pesquisa 
do Estado de Sao Paulo (FAPESP), by Conselho Nacional de e Ciencia e Tecnologia (CNPq), 
by a U.S. NSF grant INT-9805573, by Department of Energy grants DE-FG02-91ER40676 
and DE-FG02-95ER40896 and by the University of Wisconsin Research Committee with 
funds granted by the Wisconsin Alumni Research Foundation. 



14 



REFERENCES 



[1] Kamiokande Collaboration, H. S. Hirata et ai, Phys. Lett. B 205, 416 (1988); ibid. 
280, 146 (1992); Y. Fukuda et ai, ibid. 335, 237 (1994). 
1MB Collaboration, R. Becker-Szendy et al, Phys. Rev. D 46, 3720 (1992). 
MACRO Collaboration, M. Ambrosio et a/., Phys. Lett. B 478, 5 (2000); B. C. Barish, 
Nucl. Phys. B (Proc. Suppl.) 91, 141 (2001). 

Soudan-2 Collaboration, W. W. M. Allison et al, Phys. Lett. B 391, 491 (1997); Phys. 
Lett. B 449, 137 (1999); W. A. Mann, Nucl. Phys. B (Proc. Suppl.) 91, 134 (2001). 
The Super- Kamiokande Collaboration, Y. Fukuda et al, Phys. Rev. Lett. 81, 1562 
(1998); H. Sobel, Nucl. Phys. B (Proc. Suppl.) 91, 127 (2001). 
Homestake Collaboration, K. Lande et al, Astrophys. J. 496, 505 (1998). 
Kamiokande Collaboration, Y. Fukuda et al, Phys. Rev. Lett. 77, 1683 (1996). 
GALLEX Collaboration, W. Hampel et a/., Phys. Lett. B 447, 127 (1999). 
GNO Collaboration, M. Ahmann et al, Phys. Lett. B 490, 16 (2000). 
SAGE Collaboration, J. N. Abdurashitov et al, Phys. Rev. C 60, 055801 (1999); V. N. 
Gavrin, Nucl. Phys. B (Proc. Suppl.) 91, 36 (2001). 

Y. Suzuki, Nucl. Phys. B (Proc. Suppl.) 91, 29 (2001); Super-Kamiokande Collabora- 
tion, S. Fukuda et al, |hep-ex/0103032| . 

SNO Collaboration, Q. R. Ahmad et al, Phys. Rev. Lett. 87, 071301 (2001). 
K2K Collaboration, S. H. Ahn, et al, |hep-ex/ 01 03001 . 



G. L. Fogh, E. Lisi and A. Marrone, |hep-ph/0110"089i 

G. L. Fogh, E. Lisi, D. Montanino and A. Palazzo, Phys. Rev. D 64, 093007 (2001). 
J. N. Bahcall, M. C. Gonzalez- Garcia and C. Pena-Garay, JHEP 0108, 014 (2001); 
A. Bandyopadhyay, S. Choubey, S. Goswami and K. Kar, Phys. Lett. B 519, 83 (2001). 
M. ApoUonio et al, Phys. Lett. B 466, 415 (1999). 
LSND Collaboration, A. Aguilar et al, |hep-ex/010KMg . 

V. Barger, D. Marfatia and K. Whisnant, [hep-ph /0 106207| ; M. C. Gonzalez-Garcia, 
M. Maltoni and C. Pena-Garay, |hep-ph/0108073| . 

A. Piepke for the KamLAND Collaboration, Nucl. Phys. B (Proc. Suppl.) 91,99 (2001). 
V. Barger, D. Marfatia, and B. P. Wood, Phys. Lett. B 498, 53 (2001). 
The MINOS Collaboration, "Neutrino Oscillation Physics at Fermilab: The NuMI- 
MINOS Project", Fermilab Report No. NuMI-L-375 (1998). 

S. G. Wojcicki for the MINOS Collaboration, Nucl. Phys. B (Proc. Suppl.) 91,216 
(2001). 

P. G. Harris, Fermilab Report No. NuMI-L-726 (2001). 

The ICARUS Collaboration, The ICARUS Experiment, LNGS-P28/2001, LNGS-EXP 
13/89 add.1/01, ICARUS-TM/2001, see |http: / /pcnometh4.cern.ch/pubncations.htm| ; 
see also the slides from IX International Workshop on "Neutrino Telescopes", Venice, 
Italy, March (2001). 

[26] A. Rubbia for the ICARUS Collaboration, talk given at Skandinavian NeutrinO Work- 
shop (SNOW), Uppsala, Sweden, February (2001), which are available at [http:/ /pcno- 
meth4.cern.ch/publicpdf. html 



15 



[27] OPERA Collaboration, The OPERA Vj. appearance experiment in the CERN-Gran 
Sasso neutrino beam, CERN/SPSC 98-25, SPSC/M612, LNGS-LOI 8/97, October, 
1998. 

[28] OPERA Collaboration, Experiment Proposal OPERA: An appearance experiment to 

search oscillations in the CNGS neutrino beam, CERN/SPSC 2000-028, 

SPSC/P318, LNGS P25/2000, July, 2000. 
[29] V. Barger, S. Geer, R. Raja and K. Whisnant, Phys. Rev. D 63, 113011 (2001); 

B. Richter, [hep-ph/ 00082^ ; Y. Itow et al, [hep-ex/0l060T^ . 
[30] S. Geer, Phys. Rev. D 57, 6989 (1998); V. Barger, S. Geer and K. Whisnant, Phys. 

Rev. D 61, 053004 (2000). 
[31] T. Adams et a/., |hep-ph/0111030| . 

[32] D. E. Groom et al. [Particle Data Group Collaboration], Eur. Phys. J. C 15, 1 (2000). 

[33] V. Barger, K. Whisnant, S. Pakvasa and R. J. Phillips, Phys. Rev. D 22, 2718 (1980). 

[34] J. Arafune, M. Koike and J. Sato, Phys. Rev. D 56, 3093 (1997) [Erratum-ibid. D 60, 
119905 (1997)]; A. Cervera, A. Donini, M. B. Gavela, J. J. Gomez Cadenas, P. Hernan- 
dez, O. Mena and S. Rigolin, Nucl. Phys. B 579, 17 (2000) [Erratum-ibid. B 593, 731 
(2000)]; M. Freund, Phys. Rev. D 64, 053003 (2001). 

[35] A. Dziewonski and D. Anderson, Phys. Earth Planet. Inter. 25, 297 (1981). 

[36] The MINOS Collaboration, P. Adamson et al, Fermilab Report No. NuMI-L-337 (1998). 

[37] M. Messier et al, Fermilab Report No. NuMI-B-700 (2001). 

[38] H. M. Gallagher, Ph.D. thesis. University of Minnesota, 1996, available at 

[ht tp : / / hepunx. rl. ac.uk/ soudan2 /pubs/ theses. html| . 
[39] M. D. Messier, Ph.D. thesis, Boston University, 1999, available at 

|http:/ /hep.bu.edu~messier/| . 
[40] D. A. Petyt, Fermilab Report No. NuMI-L-481 (1999). We are using the new cc-efficiency 

and nc-contamination curves that were provided to us by D. A. Petyt. 
[41] M. Diwan, M. Messier, B. Viren, and L. Wai, Fermilab Report No. NuMI-L-714 (2001). 
[42] The charged current cross sections for and z/^ can be obtained in form of a table from 

|http:/ /www.cern.ch/NGS . 



[43] ICANOE Collaboration, ICANOE, A proposal for a CERN-GS long haseUne and 
atmospheric neutrino oscillation experiment, INFN/AE-99-17, CERN/SPSC 99-25, 
SPSC/P314, August, 1999, see http//pcnometh4.cern.ch/publications.html. 

[44] A. Rubbia, Nucl. Phys. B (Proc. Suppl.) 91, 223 (2001). 

[45] C. G. Arroyo, Fermilab Report No. NuMI-L-482 (1999). 

[46] V. Barger, D. Marfatia and K. Whisnant, |hep-ph/010gIMl . 



16 



FIGURES 




FIG. 1. MINOS sensitivities for the i^^ Ue and channels at the 90% C.L.. The other 

oscillation parameters are fixed at Am2i = 5 x 10~^ eV^ and sin^ 29i2 = 0.8, the best-fit solution 
to the solar neutrino problem. To determine the sensitivity to sin^ 2^13 when sin^ 2^23 = 0.85, we 
have taken 623 < 7r/4. The dashed curve would be to the left of the solid curve in the left panel if 
023 > 7r/4. 



17 



X 10 



-2 



CM 

> 



E 
<l 



0.35 



0.3 



0.25 



(V 

CD 
CN 

CM 

c 
'co 



0.9 



0.9 




sin'20 



23 



INOS 10 kt-yr 



90% C.L 
99% C.L 



theoretical input 



0.8 



10 



10 



FIG. 2. Expected allowed regions for MINOS at the 90% and 99% C.L. using Amis = 3 x 10"^ 
eV^, Amii = 5 X 10"^ eV^, sin^ 26*23 = 1, sin^ 26*12 = 0.8 and sin^ 26*13 = 0.05 as the theoretical 
input for which data was simulated. The dashed line is the Super-Kamiokande allowed region at 
the 99% C.L.. Here and in other figures, the best fit point is very close to the theoretical input 
and is consequently not shown. 
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FIG. 4. The same as Fig. | but for ICARUS. 
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FIG. 6. The same as Fig. | but for OPERA. 
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FIG. 7. The global sensitivity of the Un 



mode at 90% C.L. is shown on the left. The 



sensitivity to sin 2^23 and Am32 by combining the MINOS sensitivity in the disappearance channel 
and the ICARUS and OPERA sensitivities in the — > Vt channels at the 90% C.L. is on the right. 
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FIG. 8. Expected allowed regions from the cumulative simulated data of the three exper- 
iments at the 90%, 95% and 99% C.L. using Am^g = 3 x lO'^ eV^, Am|i = 5 x IQ-^ eV^ 
sin^ 2^23 = 1, sin^ 2^12 = 0.8 and sin^ 29i3 = 0.05 as the theoretical input. The dashed line is the 
Super-Kamiokande allowed region at the 99% C.L.. 



24 



